Chaos from turbulence: stochastic-chaotic equilibrium in turbulent convection at high 

Rayleigh numbers 
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It is shown that correlation function of the mean wind velocity generated by a turbulent thermal 
convection (Rayleigh number Ra ~ 10 11 ) exhibits exponential decay with a very long correlation 
time, while corresponding largest Lyapunov exponent is certainly positive. These results together 
with the reconstructed phase portrait indicate presence of chaotic component in the examined 
mean wind. Telegraph approximation is also used to study relative contribution of the chaotic and 
stochastic components to the mean wind fluctuations and an equilibrium between these components 
has been studied in detail. 

PACS numbers: 47.55.pb, 47.52.+j 



Introduction. Turbulent thermal convection can gen- 
erate large-scale (coherent) circulations, also known as 
mean winds. In recent decade a vigorous investigation 
of statistical properties of the thermal mean winds has 
been launched [3,0] (see for a recent review 1). The 
winds' dynamics turned out to be very complex and 
many surprising features were discovered in laboratory 
experiments and numerical simulations. Simple stochas- 
tic models @,|H were replaced by more sophisticated 
stochastic models 0,(1], which also addressed signif- 
icant three-dimensional nature of the phenomenon. In 
these models, interaction between the large-scale wind 
and the small-scale turbulence provides a phenomeno- 
logical stochastic driving term. Also some determinis- 
tic models with chaotic solutions were recently suggested 
0,01, and the idea that a large-scale instability in the de- 
veloped turbulent convection (caused by a redistribution 
of the turbulent heat flux) can be an origin of the large- 
scale coherent structures received certain experimental 
support (see Ref. [§] and references therein). 

It is a difficult problem to distinguish between stochas- 
tic and chaotic processes, especially having a developed 
turbulence as a small-scale background. Observation of 
the Lorentzian spectra for the mean winds (reported in 
the Ref. 0) makes this problem even more difficult. It 
is shown recently (see, for instance, Ref. [l(|) that many 
properties of the chaotic processes with exponentially 
decaying correlation function (i.e. with the Lorentzian 
spectra) can be reproduced by some classical models of 
stochastic processes (and vise versa). Therefore, for such 
chaotic systems stochastic phenomenological models can 
be rather useful (and relevant). However, it is necessary 
to reveal an underlying chaotic (deterministic) nature 
of the phenomenon in order to understand real cause of 
the long-term correlations in such systems. 

Chaotic mean wind. Figure 1 shows a correlation func- 
tion of a mean wind velocity, v(t), measured in a typi- 
cal thermal convection laboratory experiment in a closed 
cylindrical container of aspect ratio 1 for Rayleigh num- 
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FIG. 1: Autocorrelation function of a mean wind velocity for 
a Rayleigh-Benard convection laboratory experiment 0. The 
dashed straight line indicates the exponential decay Eq. (I). 
The insert shows a small-time-scales part of the correlation 
function defect in ln-ln scales. The straight line indicates the 
Kolmogorov's '2/3' power law for structure function. 



ber Ra ~ 10 11 (see Ref. ■ flU for description of the ex- 
periment and of the other properties of the mean wind). 
The dashed straight line is drawn in the figure in order 
to indicate (in the semi-log scales) the exponential decay 



C(t) ~ exp-(r/r ) 



(1) 



The correlation time tq ~ 600s is very large in compar- 
ison with the mean wind circulation period T c ~ 30s. 
The circulation period (turnover time) provides a time 
scale relevant to the phenomenon physics. The correla- 
tion time normalized by this time scale tq/T c ~ 20. 

The insert in Fig. 1 shows a small-times part of the 
correlation function defect. The ln-ln scales have been 
used in this figure in order to show a power law (the 
dashed straight line) for structure function: ((v(t + t) — 
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FIG. 2: The pertaining average maximal Lyapunov exponent 
at the pertaining time, calculated for the same data as those 
used for calculation of the correlation function (Fig. 1). The 
dashed straight line indicates convergence to a positive value. 



trary close trajectories of the system will diverge apart 
exponentially, that is the hallmark of chaos. At present 
time there is no theory relating \ max to Tq" 1 for chaotic 
systems. There is only a tentative suggestion that their 
values should be of the same order. To calculate \ m ax 
we used a direct algorithm developed by Wolf et al [l3j ]. 
Figure 2 shows the pertaining average maximal Lyapunov 
exponent at the pertaining time, calculated for the same 
data as those used for calculation of the correlation func- 
tion (Fig. 1). The largest Lyapunov exponent converges 
very well to a positive value X m ax — 0. 0023s -1 . 

The ambivalent nature of the turbulence-induced co- 
herent dynamics one can also see in Figure 3. This fig- 
ure shows a phase portrait reconstructed from the noise- 
reduced time series for the mean wind velocity. 

Thus, we can conclude that the thermal wind under 
study exhibit chaotic features and the long-term expo- 
nential correlation (Fig. 1) is presumably related to this 
chaotic behavior. 
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FIG. 3: Phase portrait reconstructed from noise-reduced time 
series. 



v(t)) 2 ) for r < T c 



1 - C(t) oc ((v(t + r) - v(t)) 2 ) cx r 2 / 3 



(2) 



This power law: '2/3', for structure function (by virtue 
of the Taylor hypothesis transforming the time scaling 
into the space one [13]) is known for fully developed tur- 
bulence as Kolmogorov's power law (cf Refs. H,[H). 

In order to determine the presence of a deterministic 
chaos in the time series corresponding to the velocity of 
the mean wind for the time scales T c < t < To, we calcu- 
lated the largest Lyapunov exponent : X max - A strong 
indicator for the presence of chaos in the examined time 
series is condition X rnax > 0. If this is the case, then we 
have so-called exponential instability. Namely, two arbi- 



Telegraph approximation. It is shown in Ref. 14 1 
that simple telegraph approximation of stochastic sig- 
nals can reproduce main statistical properties of these 
signals. The telegraph approximation of signal v(t) can 
be constructed as following: 



u(t) = v(t)/\v(t)\ 



(3) 



From the definition the telegraph approximation can take 
only two values: 1 and -1. Figure 4 shows a compar- 
ison between correlation function of the full signal for 
the mean wind velocity v(t) (cf Fig. 1) and correlation 
function of its telegraph approximation u(t). One can 
see very good correspondence between these correlation 
functions. Therefore, certain main statistical properties 
of the full signal can be studied using the telegraph ap- 
proximation in this case as well. 

Very significant characteristic of the telegraph signal 
is duration, r, of the continuous intervals (boxes) where 
the signal takes value 1 or value -1. Probability density 
function of the duration (or life) times r: p(t), for the 
telegraph approximation of the mean wind velocity signal 
was studied in detail in Ref. Q. It has a peak at r ~ T c . 
In the the range T c < t < tq the probability density 
exhibits a power law 



p{r) 



(4) 



with the value of the exponent a close to —1. For r > 
To the probability density p(r) decays exponentially, as 
it should be for the random telegraph signal. Figure 5 
shows cumulative probability 



P(t) = f p{r)dr 
Jo 



(5) 



In the semi-logarithmical scales the straight line indicates 
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FIG. 4: Autocorrelation functions of a mean wind velocity for 
a Rayleigh-Benard convection laboratory experiment || (solid 
curve) and of its telegraph approximation (dashed curve). 
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FIG. 5: Cumulative probability P(t) (Eq. (5)) vs. \n(t/T c ). 
The straight line indicates correspondence to Eq. (6). The 
insert shows a linear dependence of P(t) on t/T c for t < T c . 

just the power law for p(r) (Eq. (4) with a ~ —1): 

P(t) « P(T C ) + c f T- x dT « P(T C ) + c ln(t/T c ) (6) 

where P(T C ) is a constant part provided by the turbulent 
component to the cumulative probability P(t) for tq > 
t > T c . On the whole, 

P(oo) ~ AP + c ln(r /T c ) (8) 

where 

/•OC 

AP = P(T C ) + / p(r) dr (9) 

i.e. AP is contribution of the stochastic components: 
turbulence and the large-scale random telegraph signal, 



to the total cumulative probability. In our case there is a 
probabilistic equipartition of the stochastic and chaotic 
components to the total cumulative probability P(oo) = 
1: AP ~ 0.5 ± 0.04. If this probabilistic equipartition 
is universal and the constant c in the power law Eq. (4) 
is also universal (at least asymptotically), then 

1 - AP 

In r /T c = a 0.5/c (10) 

c 

is universal as well. From the present data we can esti- 
mate: c a 0.16±0.01, and ln(To/T c ) w 3 (see also below). 

Stochastic- chaotic equilibrium. This kind of probabilis- 
tic equipartition is usually related to statistical restora- 
tion of a symmetry. For instance, the spontaneous ap- 
pearing of the mean wind results in a spontaneous break- 
ing of mirror (or parity - P) symmetry. Statistical 
restoration of the parity (P-symmetry) means probabilis- 
tic equipartition of the u{t) = 1 and u(t) = — 1 events 
(Eq. (3)). This probabilistic equipartition indeed takes 
place in present case and it is clear evidence of the sta- 
tistical restoration of the P-symmetry (the directions of 
the mean wind rotation: clockwise and anticlockwise, are 
statistically equivalent). While the P-symmetry arises 
from space inversion, the T-symmetry arises from time 
reversal (reversibility). In the present case both P- and 
T-transformations result in the same mean wind switch- 
ing: u — > —u. It can be just statistical restoration of 
the combined PT-symmetry (statistical invariance under 
joint action of parity and time reversal) that results in 
the equal probability of the reversible (chaotic) and irre- 
versible (random) switchings. The statistical restoration 
of the symmetry in a mixed stochastic-chaotic motion 
indicates an equilibrium reached between the stochastic 
and chaotic components (SC-equilibrium). 

For the closed orbits in the chaotic systems the Bowen's 
theorem is an analogue of ergodic theorem. This the- 
orem provides a basis for the equality of the average over 
the phase space and the average over large-period or- 
bits. Let us now, in the terms of the Bowen's theorem 
(cf. also (HI), consider a phase volume subregion Af2: 
fi c < Ail <c (where fi c is phase volume corresponding 
to the upper turbulent scale T c ) and all the periodic tra- 
jectories passing through it. Different closed orbits can 
be distinguished by their period T. Their distribution 
does not depend significantly on location and boundaries 
of the subregion Af2. Moreover, if we consider the peri- 
odic orbits with periods in the interval (T, T + t) passing 
through the region, then we can use an estimate 



(the coarse graining is defined by the partition of phase 
space by the cells J7 C ). Let us consider an entropy cor- 
responding to the subregion AQ (cf. Ref. [lj|) and use 
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estimate Eq. (11) 



S = (3 In 
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f3 In 



(12) 



where /3 is a constant (see below, and cf. also Ref. [T7[ for 
quantum chaos) . The logarithmic growth of entropy can 
imply certain self-similarity of joint probabilistic proper- 
ties of the time series v(t) <-> u(t) (see Eq. (3)). Indeed, 
let us consider probability distribution function, p(v,t), 
for the signal v at the impulses (boxes) of length r < t of 
its telegraph approximation u. For the stationary signal 
the self-similarity has a standard form [lH ] 



p(v,t) 



(13) 



where / is certain function of the argument v/t@ (the 
statistical restoration of the parity has been taken into 
account here). Let us consider an entropy 



S(t) = dv p(v,t) lnp(v,t) 



(14) 



Changing the integration variable: v — > v/t@ and using 
equilibrium condition: S(T C ) — 0, one obtains from Eqs. 
(13) and (14): 



S(t) = p \n(t/T c ) 



(15) 



(cfEq. (12)). 

The insert in Fig. 5 shows a linear dependence of P(t) 
on t/T c for the turbulent regime that can be related to 
an analytic dependence of P(t) on t/T c for < t/T c < 
1 (the linear dependence represents the first two terms 
approximation of the Taylor expansion; the small 'jump' 
at the point t = is presumably related to transition 
from laminar to turbulent motion). In a vicinity of the 
SC-equilibrium P(t) is an analytic function of the entropy 
S(t) = /3ln(t/T c ). Applying again the Taylor expansion 
we obtain 

P{t) = P(T c ) + c S(t) + ...^P{T c )+c [3 \n(t/T c ) (16) 

where c D = (dP/dS)\ s =o (cf. Fig. 5 and Refs. 0,0 
for the transformation t — > hit). 

Discussion. It is possible that the above consideration 
is also applicable for the large-scale coherent structures 
observed in other turbulent flows (in turbulent bound- 
ary layers, for instance [1(3]). If, for instance, at the 
SC-equilibrium the exponent (3 is a global invariant 
of motion and its value is the same for t < T c and 
for t > T c , then, taking into account that the dimen- 
sionless entropy S(tq) = 1 one obtains from Eq. (15): 



ln(ro/T c ) — 1/(3 . For the Brownian processes (3 = 1/2, 
whereas for the Kolmogorov's scaling: Eq. (2), one has 
(3 = 1/3. For Kolmogorov's turbulence as a background 
it results in estimate: ln(ro/T' c ) = 3 (cf. above). If the 
background turbulence produces scaling different from 
the Kolmogorov's one, then value of the /^-exponent can 
be different from 1/3. For the Kraichnan's scaling, for 
instance [2l|, /5 = 1/4 and, consequently, ln(r /T c ) = 4 
(i.e. in the case of the Kraichnan's background turbu- 
lence the correlations related to the large-scale coherent 
structures are even more long-term than those generated 
by the Kolmogorov's background turbulence). 



The author is grateful to J.J. Niemela, to K.R. Sreeni- 
vasan for sharing their data and discussions. 
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